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Determination of (to Properties of a Solution or Solid Using Raman Ratios 

This applicat.cn .s a Linuation in part of US Patent Appl.cation Serial No. 08/814,519 

filed on March 10, 199?. 

Field of the Invention 

This invention relates to anlpparatus and a method for analyzing properties of a solution or 
solid using Raman spectroscopy and in particular to the appl.cation of Raman peak intensity 
ratios for analyzing and predicting the properties of a solution or solid. 

Background of the Invention^ 

In many industrial chemical processes, the amount of reactants, or input components, that 
are used is less than or more than the amount necessary to carry the reaction to the point 
of obtaining a desired characteristic^) of the product stream. If too little of the input 
component is used, often the desired\arget value of a characteristic from the process is 
not obtained. Alternatively, if an excessive amount of an input component is used, the 
desired characteristic may be obtained, but the excess input component is typically 
released as waste in the effluent of the process. In other cases, excessive amounts of an 
input component may cause undesirable reactions to occur that produce unwanted 
characteristics. Further, the wasted input component is economically costly and can 
become an environmental pollutant if it is released into the environment without being 
removed or recycled from the effluent. The difficulty in controlling chemical processes, 
such as bleaching, in the pulp and paper manufacturing industry can be caused by a 
number of factors including qualitative and quantitative variability of the pulp or wood 
furnish, the composition of the process chemicals, '^nd the consistency (% wood or pulp) 
of the furnish. Further, changing market requirements for paper products may require a 
paper manufacturing operation to produce a wide variety of paper grades. New paper 

\ 



1 



Patent 

51-6USCIP 

processing methods, equipment, and chemicals force the paper bleaching operation to 
adapt to these technical changes while still monitoring various characteristics of the pulp. 

It is therefore desirable to be able to precisely control the input components to obtain the 
desired target characteristic(s) with little waste. To obtain this control, a characteristic of 
the effluent of an industrial process should be precisely monitored in real time in order to 
provide feedback control on the amount of input components which should be added to 
the reactor to avoid under use or excessive use, and waste, of the input component. 

For example, in the pulp and paper industry, hydrogen peroxide and the hydroperoxy 
anion (HO/) are important input components for the oxidation and bleaching of wood 
pulps. In a typical pulp bleaching plant situation, the control of the bleaching chemicals 
is based on the brightness of the incoming pulp, the pulp flow, and the target brightness 
and pulp physical properties that are to be achieved. The factors of incoming pulp 
brightness, pulp flow, and target brightness are then used to calculate the amount of 
bleaching chemicals required to be added to the pulp to achieve a certain final target 
brightness. In another system, the brightness of the pulp is measured after bleaching 
chemicals are added and after allowing the reaction to occur for a defined reaction time. 
The resultant brightness value of the reaction is then measured and is used for feedback 
regulation of the bleaching chemicals. 

Typically with these feedback systems, the amount of hydrogen peroxide that is used 
exceeds or overshoots the amount necessary to reach a final target characteristic, such as 
pulp final target brightness, yellowness, residual peroxide, brightness efficiency, 
yellowness efficiency, and delignification efficiency. The resultant unwanted variation in 
these pulp characteristics may cause additional processing problems in the pulp and paper 
processing mill. Further, in the case of peroxide bleaching, excessive use of hydrogen 
peroxide results in waste hydrogen peroxide in the pulp effluent, which is both costly and 
environmentally harmful. 
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In order to solve these problems the prior art has offered various solutions. For example, 
United States Patent No. 4,878,998 teaches a method for bleaching of mechanical, 
thermomechanical and chemi-mechanical pulps whereby peroxide bleaching is controlled 
by the addition of a preset amount of bleaching chemicals at a first bleaching stage, 
measuring the brightness of the pulp, feed forwardly adjusting the amount of bleaching 
chemicals to be added at a second bleaching stage as a function of the measured 
brightness of the pulp from the first stage, and then bleaching the pulp at the second 
stage. 

Canadian Patent No. 2,081 ,907 teaches a method and apparatus for determining 
information characteristics of the concentration of each of at least three intermixed 
components in kraft liquors having the steps of: identifying detectable characteristics that 
are detectable in relation to the concentration of the components, developing a 
mathematical relationship between the component and the characteristics, such as 
regression analysis, analysing a sample of solution with a UV detector, and then 
controlling the concentration of each of the three components by using the information 
from the analysis of the sample. 

While current brightness sensors are able to provide a measure of the pulp brightness, 
they are unable to measure the bleaching efficiency of the bleaching reaction itself. 
Bleaching efficiency is a change in brightness of a pulp divided by the amount of 
peroxide consumed during the bleaching reactions. Further, measurement of yellowness 
efficiency, which is defined as a change in pulp yellowness divided by the peroxide 
consumed during the bleaching reactions, also requires a method by which the residual 
peroxide in the pulp effluent can be measured. Other efficiency measures relating the 
relative improvement in pulp optical and strength properties to the consumption of 
chemicals and the generation of dissolved species derived from wood are of interest, but 
are not readily available. 
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United States Patent No. 5,842,150 and WO 96/122183 to Renberg et al. disclose a 
method based on UV/VIS/NIR/IR, for the qualitative and quantitative determination of 
quality parameters in pulp and paper and/or the organic content in effluents from pulp and 
paper production by applying chemometric methods. Renberg e, al. provide a review and 

5 discussion of the "state of the art" and the need for on-line measurement of variables 
related to pulp and effluent quality. A discussion is provided for the use of spectral 
analysis with standard chemometric procedures to derive calibrated values for pulp and 
paper strength and brightness parameters and also measures of amounts of organic 
substances, for example the Total Organic Carbon (TOC), the Chemical Oxygen Demand 

1 0 (COD), and the Biological Oxygen Demand (BOD). The abstract and the disclosure 
indicate that the invention relates to UV/VIS/NIR spectroscopy, including Raman 
spectroscopy. However, the disclosure does not provide any examples or further 
disclosure with respect to the use of Raman spectroscopy. The examples disclosed in 
column 7 in United States Patent No. 5,842, 1 50 are all examples related to UV absorption 

1 5 techniques between 200 and 360 ran to some other property. It is noted that United States 
Patent No. 5,842,150 does not disclose any wavelength region outside the UV region. 
Furthermore, Raman spectroscopy is an emission technique and does not extend to 
absorption, transmittance or reflectance techniques as discussed in United States Patent 
No. 5,842,150. The reflectance technique disclosed therein is not the same as an 

20 emission by inelastic scattering as it occurs in Raman spectroscopy. The prior art does 
not relate the spectral parameters to organic indicators and does not discuss the properties 
related to the oxidative capacity of inorganic components that may exist in multiple 
oxidation states, the development of substances that contribute to scale deposition of 
effluent components, the physical properties of polymerizable species, such as the 

25 number of endgroups, the extent of network formation, and the chain length, or the 
development of bulk, yield or fiber flexibility. 

The following two references, including references within, provide a general review of 
the application and interpretation of Raman spectroscopy with respect to lignocellulosics, 
30 Chapter 9 "An Overview of Raman Spectroscopy as Applied to Lignocellulosic 
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Materials" by Umesh P. Agarwal from a book enfitled "Advances in Lignoce.lulosics 
Characterization", edited by Dimitris S. Argyropoulos, published by Ta PP i Press 1999 
(ISBN 0-89852-357-5), and an article by Umesh P. Agarwal and Sally A. Ralph m 
Applied Spectroscopy Volume 51, Number 11, 1997. PP 1648-1655, entitled "FT Raman 
Spectroscopy of Wood: Identifying Contributions of Lignin and Carbohydrate Polymers 
in the Spectrum of Black Spruce (Picea Mariana). 

Unfortunately, it is well known that there is a present lack of an appropriate method or 
device for the monitoring and control of pulp bleaching react.cn characteristics, including 
pulp final target brightness, yellowness, residual peroxide, and efficiency, with respect to 
brightness or physical strength development, the amount of peroxide consumed, or the 
loss of lignin or carbohydrate substances. Also, it is known that pH measurement probes 
and electrochemical methods of measuring hydrogen peroxide, such as the Kajaani 
Polarox sensor made by Valmet Automation, can be unreliable under pH and chemical 
1 5 concentration conditions which are typically used for pulp brightening reactions. BTG 
Spectris (Sweden) has an instrument and method of measuring the peroxide concentration 
that employs the use of a catalyst to decompose the hydrogen peroxide to generate 
oxygen gas that increases the reaction vessel pressure. This instrument, the RPA-5000, 
then relates the change in the pressure of the reaction vessel to the concentration of 
20 peroxide. This method, while providing a badly needed measure of the peroxide 

concentration, is complicated and indirect and subject to variability related to sample 
preparation and instrument maintenance. 

Measurement of the concentrations of other pulp bleaching chemicals such as sodium 
25 hydrosulfite (dithionite), chlorine dioxide, and hypochlorite present similar difficulties to 
those encountered to for hydrogen peroxide. In general these bleaching compounds are 
oxidative or reductive substances, normally existing as one or more species of inorganic 
oxianions in solution. A review of the current state of the art of pulp bleaching practices, 
including methods for measurement and control, has been published (Pulp Bleaching. 
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Principles and Practice. Carlton W. Dence, and Duglas W. Reeve Editors. Tappi Press, 
Atlanta 1996). 
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Summary of the Invention 

In accordance with the invention there is provided a method for measuring an amount of 
peroxide or peroxyl ion of a sample comprising the following steps: (a) irradiating at least 
a portion of the sample with a laser light for generating a Raman spectrum of the sample; 
(b) obtain a Raman spectrum for obtaining at least two measurements at two different 
wavenumbers, a first measurement related to a Raman intensity related to an amount of 
peroxide or an amount of peroxyl ion, the second measurement related to the other of the 
? amount of hydrogen peroxide and the amount of peroxyl ion; and (c) formulating a 

relationship between a Raman intensity for hydrogen peroxide and a Raman intensity for 
1 5 the peroxyl ion by comparing information related to the two measurements for determing 
the amount of peroxide or peroxyl ion. 



In accordance with another embodiment of the invention there is provided a method for 
determining a property of a sample comprising the steps of: (a) irradiating at least a 
3 20 portion of the sample with a laser light for generating a Raman emitted light from the 

sample; (b) obtaining at least two measurements of the Raman emitted light between 200 
cm" 1 and 4000 cm"', a first measurement at a first wavenumber and a second measurement 
at a second wavenumber; and (c) determining a non-linear relationship between the at 
least two measurements and the property of the sample. 



25 



Furthermore, in accordance with yet another embodiment of the present invention there is 
provided a method for determining a potential of an oxidative reductive process 
comprising the following steps: (a) irradiating at least a portion of the sample with a laser 
light for generating a Raman emitted light from the sample; (b) obtaining at least two 
30 measurements of the Raman emitted light between 200 cm" 1 and 4000 cm"', a first 
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1 and peroxy, ion, and a ra.io of**, reia.ed » .he m.ensi.y peaks of.be peroxide and 

^ peroxyl ion. 

1 20 I„ accordance with a further embodiment of the present invention there is provided an 
5 appar atus for determining a property of a sample comprising: a .aser hght source for 

irradiating at least a portion of the sample for generating a Raman emitted Hght from the 
sample- a detector for detecting the Raman emitted light from the sample, said detector 
for obtaining at least two measurements of he Raman emitted light, a first measurement at 
os a first wavenumber and a second measurement at a second wavenumber; and a processor 
for receiving and processing data from the detector for determinmg a non-lmear 
relationship between the at least two measurements and the property of the sample. 

,„ accordance with another aspect of the invention there is provided a system for 
30 determining a property of a sample comprising: means for determmmg a non-lmear 
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relationship between at least two measurements and the property of the sample, the at 
teas, two measurements corresponding to Raman emitted light between 200 cm"' and 
4000 cm", and the at least two measurements comprising a first measurement at a first 
wavenumber and a second measurement at a second wavenumber. 

Further more, in accordance with the invention there is provided a system for determining 
a property of a sample comprising: means for comparing at least two measurements 
including a first measurement at a first wavenumber and a second measurement at a 
second wavenumber, the at least two measurements corresponding to Raman emitted 
10 light between 200 cm" 1 and 4000 cm"' when the sample is irradiated with a laser; means 
for determining a non-linear relationship between the at least two measurements and the 
property of the sample; and, means for determining the property of the sample in 
dependence upon the non-linear relationship. 

1 5 In accordance with a further aspect of the invention there is provided a system for 
determining an amount of at least one of hydrogen peroxide and HOO- in a solution, 
comprising: means for receiving information containing data related to at least one of a 
Raman intensity peak corresponding to peroxide and a Raman intensity peak 
corresponding to peroxyl ion; and, means for processing the information to determine 

20 indicia of a concentration of at least one of peroxide and peroxyl ion, the processing 

including an analysis of at least one of data related to the intensity peak corresponding to 
peroxide, data related the intensity peak corresponding to peroxyl ion, a sum of data 
related to the intensity peaks of the peroxide and peroxyl ion, a product of data related to 
the intensity peaks of the peroxide and peroxyl ion, and a ratio of data related to the 

25 • intensity peaks of the peroxide and peroxyl ion. 
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Brief Description of the Drawings 

Exemplary embodiments of the invention will now be described in accordance with the 
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drawings in which: 

Figure 1 presents a scheme diagram of a pulp bleaching process in accordance with an 
embodiment of the present invention; 

Figure 2 shows a matrix plot for an example of a bleaching process and shows five 
dependent properties and six corresponding UV absorbance values; 

Figure 3 shows a simi.ar matrix plot, which was formed using the same charactenstic 
properties and a representative set of ratios from the absorbance values; 

Figure 4 presents an example from pulp and paper process waters and shows a plurality 
of absorbance spectra obtained at different temperatures between 20 and 80 degrees 
Celsius; 

Figure 5 shows a matrix plot with one dependent property, viz. the temperature, and 6 
corresponding UV-visible absorbance values; 

Figure 6 shows a matrix plot with one dependent property, viz. the temperature, and a 
20 representative set of ratios from the absorbance values; 

Figure 7 shows a plot presenting Raman spectra from a dispersive system with a 514.5 
nm laser and a FT system with a 1064 nm laser; 

25 Figure 8 shows a Fourier Transform (FT) scan of a Raman scattering signal from a pulp 
bleaching mixture of hydrogen peroxide, silicate, and sulfate; 

Figure 9 shows a series of Raman spectra of pressates from peroxide bleaching of pulp; 
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Figure 10 shows a matrix plot for variables related to Aspen TMP pulp bleaching with 
hydrogen peroxide, Raman intensities, Raman intensity ratios and the pulp and bleaching 
pressate properties; 

5 Figure 1 1 presents Raman Spectra of pressates from hydrogen peroxide bleaching of aspen 
TMP pulp at different pH values; 

Figure 12 shows a graph for the prediction of pulp brightness from a model based on a 
combination of Raman peak intensities and Raman peak intensity ratios; 

10 

Figure 13 resents the Raman spectra showing HOCr and HOOH peaks at 850 cm 1 and 877 
cm" 1 , respectively; 

Figure 14 shows a matrix plot of Raman intensities and Raman ratios as a function of pH; 

Figures 15 and 16 show Raman spectra of different sulfur oxianion in a dibasic form; 

Figure 17 presents spectra showing Raman intensities of a solution of sodium hydrosulfite 
(Na 2 S 2 0 4 ) oxidizing to sulfate and sulfite ions; 

Figure 18 shows the S-0 stretching region of the Raman spectrum during the oxidation of 
hydrosulfite to sulfate; 

Figure 19 presents a matrix plot showing Raman intensities, and Raman intensity ratios with 
25 time and oxidation reduction potential (ORP); 

Figure 20 presents a series of silicate Raman spectra as a function of varying hydroxide 
concentration as taken from Prabir K. Dutta and Dan-Chung Shieh, published in Applied 
Spectroscopy, Vol. 39, No. 2, pp. 343-346 (1985); 

30 
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Figure 21 shows the intensities of the different peaks as derived from the spectra presented 
in Figure 20; 

Fl gure 22 showing a matrix plot of Raman intensity ratios as a function of the HOVSi ratio; 
Figure 23 shows Raman spectra of white water samples from ANC and MWP; and 

Figure 24 shows Raman spectra of Acetic acid and Acetic acetate buffer solution. Top line: 
0.05M acetic acetate buffer solution. Bottom line: 5% acetic acid solution (scaling X 0.2). 

Detailed Description of the Invention 

Figure 1 presents a schematic diagram of a pulp bleach.ng process in accordance with an 
embodiment of the present invention. As shown in Figure 1, wood pulp 2 is bleached in a 
bleaching chamber 4, by bleaching agent hydrogen peroxide (HA), as input component. 
H 2 0 2 flows from reservoir 18 through control valve 16 into the bleaching chamber 4. 
The output of the bleaching process is a process effluent 6, which includes both the 
bleached pulp as well as the bleaching liquor. 

A sample of the effluent 6 is diverted to a Raman spectrometer 10 for the purpose of 
obtaining at least two measurements of UV light absorption of the effluent 6 and then 
calculating the ratio of the two measurements. The ratio of the Raman emitted light 
intensity measurements have been found to correlate to various characteristics of the pulp 
effluent. For example, characteristics such as pulp brightness, pH, and pulp yellowness, 
and residual peroxide can be determined through the use of different Raman 
wavenumbers in the ratio. 

The bleaching process causes structural changes in the lignin or extractive components of 
30 the pulp including ionization of the phenolic groups of the lignin molecule. Lignin 
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degradation typically results in an increase in the number of phenolic groups that can be 
ionized by changes in pH. As a result, the relative amount of ionized phenolic groups 
usually depends on the extent of lignin degradation and pH. 

In accordance with an embodiment of the present invention, the effluent 6 is filtered by a 
0.05 urn cross-flow membrane filter (Koch Filtration #5-HMF-451 PNE-PP) filter to 
remove colloidal material from the pulp effluent. If the colloidal material is not removed 
from the effluent, the resultant turbidity may interfere and offset the Raman intensities 
and hence affect the Raman ratio of the present invention. 



Computer 12 calculates a mathematical relationship between two or more Raman 
measurements, which is a ratio or combination of ratios in this embodiment of the present 
invention. Moreover, computer 12 stores in its memory a predetermined value for a 
characteristic of the effluent. 

Appropriate software provides means for receiving information pertaining to the Raman 
measurements and means for processing this information. The processing of the 
information includes the comparison of measurements, the determination of non-linear 
relationships, and the determination of a property or properties of the sample. 

Computer 12 acts as a comparing means for comparing the measured H 2 0 2 value with a 
predetermined H 2 0 2 value. This comparison step determines the actual empirical value of 
the residual weight percentage of H 2 0 2 , a characteristic of the effluent being monitored. 

25 According to this embodiment, computer 1 2 determines whether to send a feedback 
signal 14 to control valve 16 to adjust the amount of H 2 0 2 being fed into the bleaching 
process taking place in bleaching chamber 4. Thus, computer 12 together with control 
valve 16 acts as a means for adjusting the amount of input component, such as H 2 0 2; 
according to the value of the characteristic of the effluent 6, H 2 0 2 concentration, which 

30 was determined by the ratio of the Raman emitted light measurements. 
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Thus as shown in the embodiment of the invention, shown in Figure 1, the disclosed 
invention permits a real time feedback control of a pulp bleaching process. The feedback 
system provides sufficient H 2 0 2 to the pulp bleaching process without producing 
excessive waste residual H 2 0 2 or pollution. 

If desired, other input components, such as NaOH, MgS0 4 , or a chelating agent, such as 
DTPA (diethylenetriamine pentaacetic acid), are added to the bleaching reaction in 
reduced or increased amounts in accordance this feedback system. 



Richardson et al. disclose in U.S. Patent No. 5,242, 602 a method for simultaneously 
analyzing the concentration of performance indicators in aqueous systems by determining 
an absorbance or emission spectrum in a wavelength range from 200 to 2500 nm and 
applying chemometric algorithms to the absorbance or emission spectrum. The 
15 concentrations of the performance indicators determined by the chemometric algorithms 
are compared to predefined ranges for the respective performance indicators. Thus, this 
technique teaches the application of chemometrical methods for multicomponent analysts 
using multiwavelength spectroscopy. Richardson et al. describe the following four basic 
algorithms to extract and analyze features in the overall absorption or emission spectra 
20 that are specific to the qualitative and quantitative contributions from the specific 
performance indicators: 

1 . Quantification. Measurement of the absorption or emission spectra on a series of 
samples with known concentrations of the performance indicators. 

25 2. Processing. The processing of raw data to reduce noise and optimize the 

ability of the chemometric techniques to compare known spectra with unknown spectra or 
to act on specific features for the spectra of a multi-components solution to permit 
analysis of multi-components solutions or to adjust for noise or drift. The following 
preprocessing steps were explicitly identified: a) noise reduction or smoothing; b) Fourier 

30 or Walsh transformations; c) first or second derivatives; and d)correction for drift. 
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3. Analysis. Analysis of the spectra using chemometric techniques. The following 
statistical methods were specifically identified: a) principal component analysis; b) 
regression analysis including multiple regression methods; and c) discriminant analysts. 

4. Comparison. Comparison of results from the analysis to actual values. This 
step may involve the use of multiple linear regression on 2-4 Principal 
components to obtain a calibration for different performance indicators. 

In the method disclosed by Richardson et al. the chemometric algorithms are applied 
directly to the absorbance or emission spectrum. However, the Analysis methods 
described in U.S. Patent No. 5,242,602 are ineffective if there is a small variation or a 
high degree of correlation in the data matrix. The present invention discloses the use of a 
ratio or multiple ratios from the at least two measurements for computing an empirical 
value of the characteristic. This is a significant improvement from the method described 
by Richardson et al. in that the present invention can extract the variation even if the 
obtained absorbance or emission spectra are similar, i.e. they show little variation. The 
method disclosed by Richardson et al. does not provide any mention in the processing 
step for the use of ratios as a means to improve the analysis step. Furthermore, there is no 
specific mention of using ratios as a means of reducing the redundancy of the data that 
may lead to analysis problems. The invention disclosed by Richardson et al. does not 
address problems related to the reliability of the estimation of a property/characteristic 
that has a small variation relative to the total variation in the data. 

The present invention discloses the use of at least one ratio and a value formulated from 
at least two Raman measurements, each at different wavenumber, to determine a 
characteristic of a process effluent. Combinations of the disclosed ratios may be further 
used to monitor and control characteristics of the process effluent. Multiple regression 
analysis, using a forward stepwise multiple regression, was conducted using these Raman 
ratios to determine the best combination of these ratios and coefficients, (i.e. a 
predetermined relationship) which best predicted the final property of the pulp (i.e. 
optical properties and other descriptors of bleaching including delignification efficiency 
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and residual peroxide). As shown in the graphs, the present invention can use this 
predetermined relationship to determine a pulp property (final brightness, delignification 
efficiency, or residual peroxide) using a minimum number of Raman ratio measurements. 

5 A more detailed description of the problems associated with the extraction of multiple 
properties from a series of similar ultraviolet-visible spectra is given below. The use of 
ratios formulated from Raman measurements provides a means for accentuating the 
difference between the obtained spectra. A description of the problems with the 
extraction of multiple properties from a series of similar Raman spectra. The method 
1 0 disclosed by Richardson et al. could not be used to separate the variables or parameters 
since they often cause the same kind of changes in the spectra or the absolute changes are 
relatively small. The present invention makes use of relative changes to control the 
number of parameters. 
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2 1 5 The analysis methods, such as Principal Component Analysis, Regression Analysis, and 
ft Discriminant Analysis, described by Richardson et al. relate linear combinations of 

O measured variables with the observable characteristic. The instant invention relates at 

□ least one ratio and a value to a characteristic of pulp or the process effluent and implicitly 

2 accommodates nonlinear relationships between the measured values and the determined 

2 20 characteristic. Strictly speaking a chemical concentration should be linear with the 

absorbance or emission at different wavelengths. One object of the present invention is 
to obtain a relationship between a set of emission values and a characteristic that is a 
complicated nonlinear function of a plurality of different chemical and physical factors, 
such as the delignification efficiency. However, there are many problems with applying 
25 statistical analysis directly to a series of similar Raman spectra that are described below. 
Ultimately, the best result is only as good as the inherent variance in the data. Analysis 
using a set of absorbance or emission ratios alone or in combination with Raman 
emission or scattering values amounts to generating a new or expanded data set that has a 
greater variance than the original data set. A new data set will generate more accurate 
30 and stable solutions if it contains arrays that relate more closely to the characteristic to be 



15 



51-6 US CIP " Patem 

fit than the original data. Analysis using a set of absorbance or emission ratios amounts 
to generating a new data set, based on functions from the initial data set, that is used with 
the chemometric techniques. 

5 The statistical methods described by Richardson et al. require that the variation of the 
measured absorbance at different wavelengths changes substantially as a function of the 
different properties that are mathematically related to the absorbance values. Thus, small 
variation in the data matrix will result in computational difficulties in multiple-linear 
regression analysis. Further, a small variation in the relationships between variables, as 
1 0 expressed in the correlation matrix, will result in additional computational difficulties in 
the principle component analysis. In this case, small relative variations may be easily 
masked by large absolute variations. The presence of small variations contributing to the 
^ response for different variables in the presence of large variations that are similar in the 

different variables results in an ill-conditioned or singular matrix of the absorbance values 
1 5 that is not amenable to the statistical analysis described by Richardson et al. as discussed 
below in more detail using an exemplary data matrix. 
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The use of ratios obtained from ultraviolet-visible measurements can condition the data 
matrix to emphasize the variation of emission (or conversely absorbance) at one 

20 wavelength to that of another. The analysis is then performed on the relative variation 
rather than on the absolute variation. If in the process to be monitored or controlled one 
substance is transformed into another, or if a small amount of one substance can have a 
substantial effect on the concentration of another, then the method of observing the 
relative changes is much more powerful than looking at absolute changes. The present 

25 invention as defined in the claims uses at least one ratio and hence provides a means to 
selectively weigh the contributions from the at least two measurements, each at a 
different wavelengths. 

The outstanding problems associated with the analysis of similar and featureless spectra 
30 are well known and a discussion of those problems is described by O. Thomas and S. 
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Gallot in Fresenius Journal for Analytical Chemistry, 1990, Volume 338, pages 234-237; 
and by B. Karlberg in the following PCT publication WO95/01 560. The prior art does 
not provide any references to solve this problem by using a ratio. Despite numerous 
attempts to directly address the computational difficulties related to regression and 
chemometrics with colinear data, there are no reports directly addressing the 
computational problem by using ratios. 

In the case of the quantification or analysis methods described by Richardson et al., the 
data matrix would consist of a series of absorption or emission measurements on n 
different samples (indexed by /) and p different wavelengths (indexed by j). For 
example: 

Data Matrix in which the i index is the index for the sample and the j index is the index 
for the wavelength 
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In the statistical techniques described by Richardson et al., the sample correlation 
between variables is measured as a function of the variance (s fi ) and covariance (s^) of the 
data. In this case i and j are two different indexes for the wavelength for the emission (or 
conversely absorbance) measurement and k is the index for n different samples. 

n-\k=\ 
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From these calculations the sample covariance matrix 
calculated. 



S and correlation matrix R may be 
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Multiple regression calculations require the calculation of the determinant of the data 
m atrix A The determinant of an n x n matrix is the sum of the products of the elements 
of a row of the matrix and their cofactor, A cofactor ,s a determinant of a matnx 
obtained by eliminating the row and column of the element. Thus each cofactor may be 
obtained by calculating the products of e.ements of a row and the determmant of a matnx 
obtained by eliminating the row and column of each element. The ca.culat.on of the 
determinant has several important consequences: 

! Least squares multiple regression methods require calculation of the determinant 
of the data matrix to obtain a linear combination of the different variables. Co.mear data, 
or a high degree of covariance, similarity or redundancy in the data matrix will lead to an 
ill-conditioned matrix that produces unstable results. 

a) If two rows or columns are nearly identical then the determinant will 

be nearly zero. 

b) Solutions to a set of linear equations, for example as addressed when 
using multiple regression to obtain correlated physical properties from 
Raman spectra, depend on the division using the determinant of the 
matrix. Colinear variables are problematic in these types of 
calculations. Therefore, the division by a very small number leads to 
instability in the calculation in cases where the determinant is near 



zero. 
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2. In Principal Component Analysis (PCA) the colinearity problem is resolved by 
forming a new set of orthogonal variables, the principal components. The principal 
components are linear combinations of variables that express the maximum variance. 
The principal components may be used in a multiple regression technique to derive 
relationships with properties related to the variables. The regression on principal 
component solves problems related to matrix singularity leading to problems with 
inversion. However, the number of principal components may be less than the number of 
variables and some valuable information may be lost in principal components that are 
considered statistically insignificant. 

a) When the regression matrix R or the covariance matrix S is singular or 
ill-conditioned problems may occur in extracting principal 
components. 

b) If the relevant underlying effects are small in comparison with some 
irrelevant ones, then they may not appear among the first few principal 
components. The resulting component selection problem is difficult as 
an arbitrary inclusion of the first n principal components may serve to 
degrade the performance of the model. 

3. Partial Least Squares (PLS) regression is a multivariate data analysis technique 
that can be used to extracts components (now called factors) that relate several response 
(Y) variables to several explanatory (X) variables. The method aims to identify the 
underlying factors, or linear combination of the X variables, which best model the Y 
dependent variables. PLS can deal efficiently with data sets where there are very many 
variables that are highly correlated and involving substantial random noise. 

a) Richardson et al. in U.S. Patent No. 5,242,602 does not describe the use of 
PLS although other descriptions of chemometrics describe the power of 
this technique in chemometric analysis. 

b) However, the partial least squares analysis is no better than the different 
linear functions that may be described from the data. If the 
property/characteristic that is to be predicted cannot be described as a 
linear function of the variables x then the regression will introduce 
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systematic error into the regression model. 



The analysis method described in accordance with the present invention amounts to 
generating a new data matrix that has the general appearance shown below. 

Using a set of ratios the data matrix A is transformed into a new data matrix that 
may be called the ratio matrix R. This matrix is based upon a predetermined set of 
wavenumbers for the numerators and denominators of the ratios. The i index is the index 
for the sample and the j index is the index for the wavelength of the numerators, the k 
index is the index for the set of wavelengths corresponding to the set of denominators that 
match the j indexed emission (or conversely absorbance) values for each sample. The 
ratio data matrix may be generally written as: 



R- 



4,1 




\p 




A Xk x 


A ^ 


A iA 


Aj 


A i,p 


A ka, 


A Ui 


a m, 


A n,\ 


A n,j 


A n,p 


A Xa n 







The method disclosed by Richardson et al. in U.S. Patent No. 5,242,602 comprises the 
direct determination of an emission or absorbance spectrum. However, the method 
disclosed by Richardson et al. would not work with the present invention. For example, 
in order to control pulp bleaching processes it is necessary to separate the variables, i.e. it 
is necessary to control a plurality of parameters. In peroxide bleaching processes, the pH 
and the amount of peroxide are the most important parameters to control. Nevertheless, 
to control those parameters it is not sufficient to look only at the levels of those two 
parameters. It is important to look at how the reaction has progressed which means it is 
necessary to extract a measure of how much peroxide is present, to extract a measure of 
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the pH, to extract a measure of how bright the pulp is, and to extract a measure of how 
yellow the pulp is. Using the method taught by Richardson et al. it would not be possible 
to extract these parameters because often these parameters cause generally the same kind 
of changes in the measured spectra or alternatively the absolute changes in some 
5 important areas of the spectrum are fairly small. In using the relative changes, i.e. ratios, 
it is possible to extract the variables. For example, in peroxide bleaching processes there 
are two main components to consider. One component is peroxide and the other 
component is lignin. For the lignin component three aspects have to be considered: i) 
How much lignin has been released from the pulp; ii) What is the intensity of the color of 
1 0 the removed lignin components; and iii) what is the extent of lignin ionization? Thus to 
optimize this reaction the peroxide and the amount of lignin removed is maximized under 
the constraint of minimizing the color of the removed lignin. Pulp consists of about 30% 
lignin which causes its yellow color and hence removing the lignin from the pulp makes 
the pulp brighter. Lignin has a relatively broad diffuse spectrum and could be monitored 
1 5 at almost any wavelength in the UV region. Since the bleaching process involves a 
jfj plurality of parameters this process cannot simply be optimized by maximizing the 

° amount of peroxide. The problem is that the competition between productive bleaching 

reactions involving the nucleophilic peroxyl anion compete with destructive reactions 
initiated by a (second order) bimolecular decomposition involving hydrogen peroxide and 
20 the peroxyl anion. Variation in pH or peroxide concentrations yields a point where the 
peroxide decomposition dominates over the productive bleaching reactions. Before that 
point the improvements in the bleaching process are basically proportional to the pH or 
the peroxide concentration until a maximum is reached in both. At the maximum the 
peroxide decomposition leads to darkening reactions that neutralize the productive 
25 bleaching reactions. Thus at this pH and this peroxide concentration the amount of 
bleaching that has been done is no longer related to the amount of peroxide consumed 
and also it is no longer related to the amount of lignin released. In fact, the trend is 
opposite and has a non-linear equation. In the beginning, the amount of bleaching is 
linear for a certain region when observing the bleaching versus the amount of peroxide 
30 consumed or the bleaching versus the amount of lignin released. The bleaching is linear 
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in some regions but in fact it is a parabolic function and one observes a maximum. Using 
the method disclosed by Richardson et al., i.e. a linear function, it appears that the 
bleaching process is progressing but in reality this is not happening. Thus, it is necessary 
to identify when an unproductive loss in hydrogen peroxide and an unproductive 

5 generation of the lignin occurs. This becomes important when the peroxide starts to 
decompose and superoxide and hydroxyl free radicals are formed. These free radicals 
start reacting with the lignin and large amounts of lignin are removed. However, the pulp 
whiteness decreases and increased amounts of peroxide are used without reaching 
optimum bleaching conditions or improving the brightness of the pulp. Therefore, it is 

10 important that this point is characterized. The present invention discloses that when this 
point is approached, the color of the lignin changes with a resulting increase in the 
absorbance value at 350 nm, for example. The relative amount of color in the lignin 
increases substantially at that critical point. The method taught by Richardson et al. 
would detect this increase in absorbance as a result of lignin removal. Lignin has 

1 5 intrinsic or natural color and when more lignin is removed more color is observed. This 
means, that in the initial bleaching stages lignin is removed with the intrinsic color. 
However, when the peroxide starts to break down, then a new kind of reaction occurs and 
more lignin is removed that is highly colored. Problems arise because at the same time 
these highly colored lignin species are generated on the pulp. The method by Richardson 

20 et al. would measure changes in the amount of lignin removed but because this method 
does not make use of relative amounts it cannot distinguish between the desired removal 
of natural lignin and the undesired case where removal of highly colored lignin occurs as 
a result of a competing reaction when the peroxide breaks down. Thus without making 
use of absorbance or emission ratios it is not possible to obtain sufficient information to 

25 distinguish between one reaction happening and another reaction happening. The method 
taught in U.S. Patent No. 5,242,602 can not manage this non-linear relationship. 

The following example demonstrates the advantage of the present invention and how it 
addresses the critical point of extracting the variation by using ratios instead of applying 
30 analysis techniques directly on the spectra or on the processed spectra as it was taught in 
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the method disclosed by Richardson et al. The presented example utilizes ratios obtained 
from UV absorbance measurements. This example is an example of a bleaching process 
and shows five dependent properties and six corresponding UV absorbance values. The 
absorbance values are representative, but are chosen to provide a large variance. Figure 
2, shows a matrix plot for an example of a bleaching process and shows five dependent 
properties and six corresponding UV absorbance values. Each miniature plot in Figure 2 
shows a bivariate plot of the corresponding diagonal elements. The variable on any given 
plot axis is determined by vertical and horizontal tracing to the variable at the plot 
diagonal, much like a mileage chart on a map. This matrix plot demonstrates that the 
absorbance functions are highly colinear with each other. There are five different non- 
linear relationships between the absorbance values and the different properties: pH, 
residual H 2 0 2 > brightness, yellowness (indicated as BSTAR) and brightness efficiency. 
Furthermore, it is seen that there is scatter in the direct relationships between some of the 
characteristic properties and the absorbance values. 

Figure 3 shows a similar matrix plot, which was formed using the same characteristic 
properties and a representative set of ratios from the absorbance values. Note that the 
colinearity between the ratio variables is much smaller than the colinearity between the 
absorbance values. Furthermore it is seen that the different UV-visible ratios vary 
linearly with different properties. There is a high correlation coefficient between 
different single ratios and each of the characteristic properties. These correlation 
coefficients are given in Table 1 below. 
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Table 1 



Property 


UV ratio 


Ratio 

Correlation 


Wavelength 
UV 

absorbance 


Single Linear 
Correlation 


pH 


A230/A250 


-0.98 


A250 


0.96 


H 2 0 2 Residual 


A230/A280 


0.96 


A280 




Pulp Brightness 


A350/A250 


-0.83 


A400 


-0.80 


Yellowness (B*) 


A350/A280 


0.95 


A400 


0.95 


Brightness 
Efficiency 


A230/A400 


0.52 


A450 


-0.42 



Principle Component Analysis: 

Using 9 UV values, only one factor may be extracted. This factor explains 97.6% of the 
data variation. There is one other factor with an eigenvalue between 1 .0 and 0. 1 . 

However, using 9 UV ratios, two factors may be extracted. One factor explains 67.6% of 
the data variation and the other explains 22.2% of the data variation. There are three 
other factors with eigenvalues between 0.1 and 1.0. 

The ability to extract more factors using principle component analysis indicates that the 
data variation may be broken down into more independent factors using the ratio method. 

Figure 4 presents an example from pulp and paper process waters. Figure 4 shows a 
plurality of absorbance spectra obtained at different temperatures between 20 and 80 
degrees Celsius. The small variation in the obtained spectra shows that there is a 
chemical change as a function of the temperature. These formless spectra with their small 
variations as a function of the process variable are typically encountered in the pulp and 
paper industry as well as in many other industries consuming a large amount of process 
water. Figures 5 and 6 present the matrix plots for this example. The set of spectra for 
this example shows small variations in the absorbance values. 
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Figure 5 shows one dependent property, viz. the temperature, and 6 corresponding In- 
visible absorbance values. In this case, one principle component accounts for 97% of the 
data variation. The general trend is for the absorbance values to decrease with increasing 
temperature. The principle components method can only resolve this effect, even though 
it is clear that more is happening. A second principle component may be found if 
eigenvalues between 0.1 and 1.0 are allowed. This second eigenvalue accounts for only 
2% of the variation. 

Figure 6 shows one dependent property, viz. the temperature, and a representative set of 
ratios from the absorbance values. In this case, two principle components account for 
61.2% of the data variation and 27.6% of the data variation. A third principle component 
may be found if eigenvalues between 0.1 and 1.0 are allowed. Once rotated, the three 
eigenvalues may account for 40%, 34%, and 25% of the data variation. 

In accordance with an embodiment of the present invention it is taught that an efficiency 
and ultimate extent of a mechanical pulp bleaching operation using hydrogen peroxide is 
monitored and controlled by analysis of the UV-visible spectrum or Raman spectrum of 
the pressate from the bleaching process. It was not previously known that the reactions 
on the components in the pressate were representative of the pulp bleaching process. 
Bleaching control has always been done using pulp brightness and without observing the 
color of the pressate. Measurements of pulp brightness suffer from variability related to 
pulp consistency, pulp surface area or scattering and insensitivity. The present invention 
teaches that the amount of lignin in the pressate and the color of the pressate, as 
determined by UV ratios, reflected the bleaching progress during a mechanical pulp 
bleaching process. Previously, these bleaching processes have been described only in 
terms of the reactions on the pulp and the reactions on the components in the pressate 
were thought to be trivial and different from the bleaching reactions. The method of 
using a series of ratios of absorbance or alternatively emission values from pressates has a 
further advantage because pressates are easier to analyze optically than the pulp. At the 
point of maximum pulp brightness the response of pulp brightness is a flat function of 
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pulp bleaching variables such as pH and peroxide concentration. For this reason, pulp 
brightness is a very poor control parameter. Another unexpected and advantageous result 
is that a combination of UV ratios or Raman ratios provides sensitive measures at this 
point where traditional bleaching control sensors that monitor a pulp brightness fail. The 
analysis disclosed in the patent describes that hydrogen peroxide bleaching of mechanical 
pulps can be controlled by monitoring a combination of the amount of hydrogen peroxide 
anion, the amount of lignin removed, and the relative amount of color removed during the 
bleaching process. Three points may be made concerning these parameters: 

1. The removal of lignin was identified as an important process in the bleaching 
of pulp. In the prior art this was not thought to be the case for mechanical pulps. 
However, excessive lignin removal can occur during inefficient bleaching processes that 
leads to a parallel increase in colored components. The best bleaching occurred when the 
most lignin could be removed while generating the smallest amount of colored 
components. 

2. The bleaching process was most efficient when the relative amount of 
hydrogen peroxide anion was maintained at high levels. Underlying this simple 
statement are the realizations that the principle cause of inefficient bleaching is the 
degradation of peroxide at high concentrations and high pH values and an important 
reason for incomplete bleaching is low pH values that do not fully activate the hydrogen 
peroxide to the hydrogen peroxide anion. Furthermore, peroxide degradation leads to 
pulp yellowing and darkening processes. The optimum bleaching process then is one at 
which the bleaching occurs at the highest levels of peroxide anion concentration that are 
consistent with the beneficial effects of bleaching but which exclude the negative effects 

of peroxide degradation. 

3. A strategy of maximizing the relative concentration of the peroxide (HOO) 
bleaching agent and the beneficial effects of the bleaching agent of color removal are 
taught in accordance with an embodiment of the present invention. This method does not 
require a multiple regression with respect to a concentration in the effluent. A learning 
step with respect to the final pulp properties is advantageous but not necessary. 
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As of yet, reliable means for determining hydrogen peroxide concentration are not 
available. The value of hydrogen peroxide consumed at a paper mill is typically between 
two and ten million dollars a year. Much of this peroxide is wasted and the development 
of a reliable means for monitoring hydrogen peroxide in the pulp and paper industry is 
5 important in making the bleaching process more economical. The advantages of using a 
set of UV absorbance or alternatively Raman scattering ratios obtained from a filtrate to 
determine the peroxide concentration have not been previously recognized. A simple 
measurement of the hydrogen peroxide ion concentration as a function of UV absorbance 
or Raman emission would be obvious. However, the measurement of the hydrogen 
1 0 peroxide concentration in a background of strongly UV absorbing substances, such as 
lignin, presents substantial deconvolution problems. Furthermore, the measurement of 
the peroxide concentration under pH conditions where the relative amounts of hydrogen 
peroxide and the hydrogen peroxide concentrations are varying is a difficult problem. 
The present invention obviates this problem, in that by using ratios, i.e. relative values, 
1 5 obtained from ultraviolet-visible measurements the contributions from lignin and the 
Ul variation due to pH are factored out. 

^ A plurality of different methods for measuring peroxides are known. Among those, 

D polarigraphic measurements are thought to be unreliable and commercial methods using 

20 catalytic peroxide decomposition were not yet available. Presently, a hydrogen peroxide 
measurement is available from BTG a Division of Spectris Technologies, but the 
reliability of this measurement is still in question and it has not been widely adapted. 

In accordance with the present invention ratios of Raman peak intensities are used to 
25 predict the properties of a solution or a solid such as pulp that is processed with the 
solution. The method and apparatus of the present invention extend the multiple ratio 
strategy to Raman scattering intensities observed at visible or near IR wavelengths. In 
this case the intensity of the Raman shifted data is used to create the ratio. These 
intensities are related to the concentration of species dissolved in the liquid. A preferred 
30 measurement is using a transmissive Raman scattering measurement using a Nd: YAG 
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(1064 nm) laser to minimize sample fluorescence. If desired, baseline corrected spectra 
and a scattering from a water reference is subtracted before extraction of intensities for 
ratios. The Raman scattering intensities provide a good measure of the concentration of 
small oxygenated molecules. Relevant small molecules and complex ions in the pulp and 
5 paper industry include, but are not exhaustive, S0 4 2 ', S0 3 2 \ H 2 0 2 , C102, HC103, silicates, 
acetic acid, Chloric Acid HC103, Chlorate C103(-) s Chlorous Acid HC102, Chlorite 
C102(-), Hypochlorous Acid HCIO. Hypochlorite CIO(-), phosphate, nitrate, nitrites. If 
desired, the present invention is used to measure large molecules such as hemicellulose, 
extractives and pectic substances. The use of Raman spectra to measure small, 
10 oxygenated molecules is well documented and extensive lists exist of Raman shifts for 
molecules and functional groups, which are found for example in "Infrared and Raman 
Spectra of Inorganic and Coordination Compounds Part A: Theory and Applications in 
Inorganic Chemistry" Kazuo Nakamoto, 1997, John Wiley and Sons, Inc, New York. 
The value of using Raman emissions to measure and control species of oxidative and 
5 1 5 reductive oxygenated inorganic species has not been recognized before this invention. 
% Furthermore, the importance of measurement of the relative concentrations of different 

0 related species, for example: HOOH and HOO- or C10 2 , C10 2 \ HCIO and CIO; or S 2 0 4 2 \ 

fi: 

0 S0 3 2 ", S0 4 2 ", has not been recognized or adapted. 



pi 

tea? 



20 Figure 7 shows an example of dispersive and Fourier Transform (FT) Raman spectra to 
3 demonstrate an interference using a laser at shorter wavelengths. The plot presents 

Raman data from a dispersive system with a 514.5 nm laser and a FT system with a 1064 
nm laser. The data from the dispersive system shows interference that probably is due to 
fluorescence, at 800-1000 cm" 1 . With advances in Raman spectroscopy lasers at a longer 
25 wavelength than 514.5 nm are used to obviate a fluorescence problem. Thus a preferred 
configuration in accordance with an embodiment of the present invention is a laser 
operating at a higher wavelength than 5 14.5 nm. 

Figure 8 shows a Fourier Transform (FT) scan of a Raman scattering signal from a pulp 
30 bleaching mixture of hydrogen peroxide, silicate, and sulfate. Steps in the processing of 
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the raw data include a fast Fourier transform, a baseline correction and then subtraction of 
the water signal. The peak at 400 cm' 1 and below is an experimental artifact due to a 
silicate cell and the detection geometry. An optimized system would have a different 
geometry, such as 90 degree detection, or cell material, such as sapphire, to minimize the 
interference of the cell material when measuring silicate solutions and colloids. 

Figure 9 shows a series of Raman spectra of pressates from peroxide bleaching of pulp. 
These spectra were obtained using a 1064 nm laser with FT signal processing, baseline 
correction and subtraction of water spectra. The samples were quantitatively diluted to 
pH 7 before their measurement. The series progresses from high bleaching pH at the top 
to low bleaching pH at the bottom. Raman peaks are observed at approximately 530 cm' 1 
for silicate, at approximately 877 cm' 1 for hydrogen peroxide (H 2 0 2 , at approximately 
990 crnl for sulfate, and at approximately 1077 cm' 1 for and carbonate. The Raman 
peaks for sulfate shown in Figure 8 appear due to the addition of sulfate. Sulfate was 
added in order to lower the pH value. The Raman peaks for carbonate are observed 
because C0 2 is absorbed from the surrounding atmosphere to form carbonate ions. 
Alternatively, this method is used for predicting scale formation since it is capable of 
detecting sulfate ions and carbonate ions. In this case, care is taken that the pressates are 
not exposed to carbon dioxide. 

Table 2 shows a relationship between pulp yield and Raman intensity ratios. 



Table 2 

Regression Summary for Dependent Variable: YIELD 
R= .99912365 R 2 = .99824807 Adjusted R 2 - .99562018 



FH 2^1=379 87 d<.00263 Std.Error of estimate: .27325 








St. Err. 




St. Err. 








BETA 


of BETA 


B 


ofB 


t(2) 


p-level 


Intercpt 






37.37 


4.841 


7.72089 


.016364 


RH202U 
V2 


.992604 


.080541 


8322.67 


675.313 


12.32417 


.006520 


H202 


.958102 


.091885 


32733.36 


3139.227 


10.42720 


,009072 


RSILH20 


.942172 


.122940 


38.20 


4.985 


7.66366 


.016604 
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Table 3 shows a relationship between residual hydrogen peroxide and Raman ratios. 



Table 3 



Regression Summary for Dependent Variable: H202 residual 
R= .99822460 R 2 = .99645235 Adjusted R 2 ~ .991 13086 
F(3 ,2)=187.25 p<.00532 Std.Error of estimate: .1 1298 







St. Err. 




St. Err. 








BETA 


of BETA 


B 


ofB 


t(2) 


p-level 


Intercpt 






-12.085 


1.7914 


-6.74606 


.021275 


RH202U 
V2 


1.045293 


.122605 


2546.622 


298.7004 


8.52567 


.013480 


RH202SI 
L 


1.376609 


.171798 


1.718 


.2144 


8.01298 


.015220 


RSILH20 
2 


1.395214 


.236563 


16.437 


2.7869 


5.89786 


.027565 



Figure 10 shows a matrix plot for Aspen TMP pulp bleaching with hydrogen peroxide 
and the pulp and bleaching pressate properties. This matrix plot shows variables 
describing the peroxide bleaching process and Raman peak intensity and intensity ratios 
from bleaching pressates obtained by bleaching with hydrogen peroxide under varied pH 
conditions. The samples were diluted to a constant pH before measurement so that the 
information in them relates to changes developed during the bleaching process. The pulp 
variables include an ultimate brightness, bulk, which is a measure of the specific volume 
of the pulp, and yield, which is the % mass of the bleached pulp relative to the 
unbleached pulp. The bleaching pressate values are characterized by the process variable 
pH, the residual peroxide as a % on pulp when 4% by weight of peroxide and pulp that 
was initially put on the pulp, and the total dissolved solids (TDS). Raman intensity ratios 
provided relate to the peak intensity for silicates at -530 cm* 1 (silicate, -530 cm' 1 ), and 
hydrogen peroxide at 877 cm' 1 (H 2 0 2 , 877 cm 1 ). The Raman intensity ratio of 
RH202SIL presented in Figure 10 is a ratio of (H 2 0 2 (877 cm" 1 ) intensity/ intensity of 
silicate (-530 cm" 1 ) and intensity of Raman intensity to UV absorbance is given by the 
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following ratio RH202UV2 as expressed by (H 2 0 2 (877cm ')/UV absorbance at 280 nm). 
For each subplot the X-axis relates to the relative concentration of the variable above and 
the Y-axis relates to the relative concentration of the variable to the right. 

The ratio formulated between the Raman peak intensity for the peroxide and the silicate is 
of importance for indicating an optimal stabilizing effect through the silicate. Silicate is 
added to the bleaching solution to stabilize the peroxide. Thus the amount of silicate with 
resoect to hydrogen peroxide is maximized such that it provides best stabilizing effects 
for the amount of hydrogen peroxide. 

Table 4 shows data for peroxide bleaching at different pH values. The Raman intensity 
ratios are compared to pulp and pressate properties. 



Table 4 



Start 
pH 



Pulp and Pr essate Properties 

rightne 



Bleach 
Yield 



residual 



g/i 



TUIp" 



ISO 



9.7 

10.0 

10.3 

10.5 

10.7 

11.1 

11.5 

12.1 



96.1 
96.0 
93.6 
92.8 
93.0 
90.2 
88.9 
86.7 



2.38 
2.82 
1.90 
1.89 
1.26 
0.91 
0.04 



ITT 
1.16 
1.76 
1.95 
2.21 
2.20 
2.69 
3.15 
5.26 



74.52 
77.97 
78.93 
78.99 
79.68 
81.82 
83.45 
80.33 



Silicate 
Bulk 530 cm-1 
T27 



Raman intensitie s of pH adjusted pressate 

877 cm-1 990 cm-1 1077 cm-1 
u~35" 



3.17 
3.12 
2.80 
2.92 
2.38 
2.60 
2.00 
1.69 



0.28 
0.32 
0.21 
0.20 
0.32 
0.30 
0.40 
0.41 



1.05 
1.08 
1.06 
1.05 
0.98 
0.98 
0.86 
0.89 



1.27 
1.29 
1.24 
1.24 
1.34 
1.28 
1.46 
1.38 



0.92 
1.04 
0.97 
0.96 
0.91 
0.94 
1.01 
0.97 



Figure 1 1 presents Raman Spectra of pressates from hydrogen peroxide bleaching of 
aspen TMP pulp at different pH values. The spectra are obtained directly, i.e. without 
dilution or pH adjustment. The Raman peaks for hydrogen peroxide (H202) at 877 cm" 1 , 
for carboxylic acid (COO-) at 924. cm' 1 , for a C-H bending mode at 1350 cm ', and for a 
C-H vending mode at 1415 cm' 1 are shown. The relative decrease of the hydrogen 
peroxide peak and increase in the peaks representing different bleaching by-products can 
be related to pulp properties developed during bleaching. During the pulp bleaching 
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process with hydrogen peroxide competing chemical processes result in the destruction of 
colored species, the cleavage and removal of colored species, and the formation of 
colored species during non-productive cleavage of lignin substances. Peroxide 
concentrations decrease during both productive and non-productive reactions through 
bimolecular degradation of the peroxide or reaction with the wood substances. The 
spectra in figure 1 1 and the associated tables 5 and 6 below, show that the productive 
reactions leading to high brightness with minimal loss of yield and bulk generate 
carboxylic groups with Raman intensities at 924 cm '. Smaller amounts of brightness 
gain and greater yield or bulk loss correlate well with the increase in the intensity of 
Raman emission peaks at 1350 cm and 1415 cm Hence the greatest improvement in 
brightness gain with the least lost of yield or bulk occurs when the amplitude of the 
intensity ratio of the peak at 924 cm to the intensity at 1 350 cm or 141 5 cm is 
maximized. The peak intensity for the hydrogen peroxide at 877 cm should be 
maximized to preserve the hydrogen peroxide and to prevent unnecessary radical 
reactions. These results provide a basis for improved relationships between the Raman 
intensities and the pulp properties shown in Tables 7-12. In Table 7 the brightness 
development is optimized when the relative consumption of the hydrogen peroxide to the 
development of the C-H peak at 1350 cm -1 is minimized. In Table 9 the relates the bulk 
to the product of the peaks at 1350 cm -1 and 1415 cm -1 as formation of compounds 
contributing to these both of these peaks contribute to brightness loss. 



Table 5 

PH ISO % 

initial final brightness 



9.3 


6.4 


72.17 


9.6 


6.5 


72.56 


10 


7.1 


74.49 


10.2 


7.3 


75.17 


10.5 


7.4 


76.05 


10.7 


7.6 


76.03 


11 


8 


77.58 


11.4 


8.4 


78.32 


11.8 


9 


78.25 



H2o2 

eness bulk % in pressate 

164.5 2.931 0.577 

149 2.967 0.548 

148 2.901 0.359 

143 2.774 0.311 

147 2.750 0.277 

152 2.815 0.245 

142 2.620 0.147 

132 2.453 0.088 

129 1.954 0.075 



% on pulp TOC 




yield 


2.974 


1580 


98.102 


2.825 


1550 


98.138 


1.850 


3030 


96.361 


1.604 


3610 


95.664 


1.428 


4070 


95.112 


1.262 


4570 


94.511 


0.760 


6550 


92.133 


0.451 


6400 


92.313 


0.388 


13100 


84.266 
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Table 6 

PH 
initial 
9.3 
9.6 
10 
10.2 
10.5 
107 
11 

11.4 
11.8 



final 



6.4 
6.5 
7.1 
7.3 
7.4 
7.6 
8 

8.4 
9 



UV ABS 
280 nm 
0.588 
0.578 
0.867 
1.013 
1.183 
1.269 
2.164 
3.193 
4.371 



Raman Intensities 



I 



_876 

1.960 
1.585 
1.207 
0.996 
0.888 
0.797 
0.656 
0.324 
0.587 



I924 I 
0.532 
0.434 
0.779 
0.893 
0.964 
0.971 
1.559 
1.647 
1.782 



1350 



I. 

0.364 
0.263 
0.334 
0.441 
0.469 
0.471 
0.731 
0.877 
1.216 



1415 
0.496 
0.409 
0.666 
0.653 
0.797 
0.795 
1.181 
1.334 
1.750 



Figure 12 shows a graph for the prediction of pulp brightness from a model based on a 
combination of Raman peak intensities and Raman peak intensity ratios. The observed 
values are plotted against the predicted values with the ISO brightness being a dependent 
variable and the independent variables being the Raman intensities and intensity ratios. 

Tables 7 to 12 below demonstrate the prediction of pulp properties from Raman ratios 
and Raman intensities. The pulp properties presented are brightness in Table 7, bulk in 
Tables 8 and 9, residual hydrogen peroxide on the pulp in Table 10, and total organic 
carbon (TOC) in Tables 1 1 and 12. 



Table 7 

Regression Summary for Dependent Variable:Brightness 
R= .99774854 R 2 = .99550215 Adjusted R 2 - .99280344 
F(3,5)=368.88 p<.00000 Std.Error of estimate: .19239 



Raman 



Intensity 



Intercept 



1876/11 35 
O 



11415 



1876 



BETA 



-.318581 



.279821 



-.444773 



St. Err. 



of BETA 



.111927 



.054338 



B 



77.16 



-.35 



.102212 



750.02 



-1016.37 



St. Err. 



ofB 



.4319 



.1226 



t(5) 



178.6687 



145.6442 



233.5693 



-2.8463 



5.1497 



-4.3515 



p-level 



.000000 



.035979 



.003616 



.007349 



Table 8 

Regression Summary for Dependent Variable: BULK 
R= .98249292 R 2 = .96529234 Adjusted R 2 = .9603341 1 
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F( 1 ,7)=1 94.68 p<.00000 Std.Error of estimate: .06333 







St. Err. 




St. Err. 








BETA 


of BETA 


B 


of B 


t(7) 


p-level 


Intercpt 






3.264 


.04657 


70.0903 


.000000 


11350 


-.982493 


.070415 


-521.783 


37.39594 


-13.9529 


.000002 



Table 9 

The Ratio Predicts bulk better 

Regression Summary for Dependent Variable: BULK 
R= .991 13735 R 2 = .98235324 Adjusted R 2 = .97983227 







St. Err. 




St. Err. 








BETA 


of BETA 


B 


ofB 


t(7) 


p-level 


Intercpt 






3. 


.022 


138.9712 


.000000 


P1350141 


-.991137 


.050209 


-127757. 


6471.954 


-19.7401 


.000000 



Ul 
□ 

5 

s« : 

\ y 



Table 10 



10 



15 



Regression Summary for Dependent Variable: H202 residual on Pulp 
R= .99806738 R 2 = .99613850 Adjusted R 2 = .99382159 



1 \J,~>J 


y ' y -.www 


St. Err. 




St. Err. 








BETA 


of BETA 


B 


ofB 


t(5) 


p-level 


Intercpt 






.799 


.16514 


4.84114 


.004710 


R876135 
O 


.360850 


.103708 


.163 


.04689 


3.47948 


.017669 


11415 


-.252501 


.050347 


-279.289 


55.68855 


-5.01519 


.004052 


1876 


.427770 


.094706 


403.384 


89.30762 


4.51680 


.006301 



Table 11 

Regression Summary for Dependent Variable: TOC 
R= .96431621 R 2 = .92990576 Adjusted R 2 = .91989230 
F( 1 ,7)=92.866 p<.00003 Std.Error of estimate: 1002.7 



Intercpt 
11350 



BETA 



.964316 



St. Err. 
of BETA 



.100067 



B 

-1394. 
5705570. 



St. Err. 
ofB 
737.4 
592068.1 



J(7) 

-1.89022 

9.63668 



p-level 

.100643 

.000027 



Table 12 



20 Regression Summary for Dependent Variable: TOC 
R= .99132396 R 2 - .98272319 Adjusted R 2 = .9723571 1 
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F(3,5)=94.802 p<.00008 Std.Error of estimate: 589.00 







St. Err. 




St. Err. 








BETA 


of BETA 


B 


ofB 


t(5) 


p-level 


Intercept 






-1989. 


1509. 


-1.31791 


.244683 


Product 
I1350«I14 

15 


2.031475 


.447598 


29173073 
E2 


64277456 
6. 


4.53862 


.006177 


Ratio 
1924/11 35 
0 


.208308 


.072085 


2406. 


832. 


2.88976 


.034197 


UV280 
abs 


-.993688 


.435777 


-2691. 


1180. 


-2.28027 


.071514 



Figure 13 presents the Raman spectra showing HOO" and HOOH peaks at 850 cm" 1 and 
877 cm" 1 , respectively. The spectra presented in Figure 13 are shown after subtraction of 
a reference spectrum of water and baseline correction. When looking at the Raman 
spectrum of hydrogen peroxide two peaks are observed for peroxides. The intensity of 
the two peaks varies with the pH value of the solution, i.e. at a lower pH value more 
hydrogen peroxide is observed and at a higher pH value more peroxyl anion (HOO") is 
observed as a result of ionization of the hydrogen peroxide. In bleaching processes it is 
important to know what percentage of the hydrogen peroxide is ionized at a particular pH 
value. This is determined by looking at the two Raman peaks at 850 cm" 1 and 877 cm" 1 
and comparing them with each other. The active species in bleaching processes are the 
ionized species and thus the amount of peroxyl anion (HOO ) needs to be maximized 
since too much peroxyl anion (HOO ) would cause a degradation of the hydrogen 
peroxide. In accordance with an embodiment of the present invention Raman 
spectroscopy conveniently provides a means for determining an extent of ionization of 
hydrogen peroxide by comparing the Raman peaks for the peroxyl anion and hydrogen 
peroxide. 

Figure 14 shows a matrix plot of Raman intensities and Raman ratios as a function of pH. 
The considered Raman intensities are the Raman intensities for hydrogen peroxide 
HOOH (1877) and peroxyl ions HOO' (1850) as a function of pH. Two ratios of Raman 
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peak intensities are also shown in the matrix plot. The Raman intensity ratio of R850877 
is the intensity of the peak at 850 cm' 1 to the intensity of the peak at 877 cm" 1 , the Raman 
intensity ratio R877T is the ration of the 877 cm"' peak to the sum of the two peak 
intensities at 850 cm' 1 and 877 cm' 1 . These intensity ratios give new non-linear variables 
that relate to the peroxide solutions properties. The matrix plot shows the relationship 
between the brightness and the Raman peak intensities of the HOO- and HOOH peroxide 
species. The use of these peak intensities for the measurement and control of hydrogen 
peroxide concentrations in pulp and paper and water treatment has not been described. 
Thus in accordance with an embodiment of the present invention the Raman peak 
intensities of the Raman peaks at 850 crn 1 and 877 cm -1 are used to control the hydrogen 
peroxide concentration. The ratio of the peak intensities provides a new relationship with 
respect to the solution properties as shown in plot 1 4C. The ratio is different than the 
direct peak intensities as shown in Figures 14F and 14H. Raman peak intensity ratios or 
products provide different means the properties are modeled by influencing variables.A 
second ratio may be defined by the intensity of one peak to the sum of the intensity of 
both Raman peaks. This ratio is expected to follow relationships relavant to the pKa of 
the peroxide. 

The matrix plot presented in Figure 14 does not show a product of the intensities of the 
two Raman peak intensities at 850 cm' 1 and 877 cm 1 . However, the product of these two 
Raman peak intensities provides similar results as the ratio between two Raman peak 
intensities. 

The Measurement of the Redox Properties and Relative Concentrations of Related 
Oxidative and Reductive Species 

Figures 15 and 16 show Raman spectra of different sulfur oxianion in a dibasic form. 

Figure 17 presents spectra showing Raman intensities of a solution of sodium 
hydrosulfite (Na 2 S 2 0 4 ) oxidizing to sulfate and sulfite ions. The management of 
oxidizing and reducing substances in industrial applications is problematic in part 
because measures such as an oxidation reduction potential (ORP) are very sensitive to 
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pH, ionic strength, temperature and the influence of interfering substances. The use of 
Raman peak intensities provides a means to directly measure the concentrations and 
relative concentrations of the different species contributing to the oxidation potential of 
the solution. 

Figure 18 shows the S-0 stretching region of the Raman spectrum during the oxidation of 
hydrosulfite to sulfate. The isobestic point is the point at which the total concentration is 
always a same function of intensity no matter what proportion of the two different 
species, i.e. hydrosulfite and sulfate, are present. 

Figure 19 presents a matrix plot showing Raman intensities, and Raman intensity ratios 
with time and oxidation reduction potential (ORP). Regression lines marked in the 
matrix plot of Figure 19 demonstrate that the Raman intensity ratios correlate better with 
the ORP than their components. Plots P and R within Figure 1 9 demonstrate that a 
Raman intensity ratio provides the best direct measure of the ORP of the sulfur species in 
solution since these plots are almost linear. 

Table 13 provides a regression summary for the oxidation reduction potential and Table 
14 shows the correlations for hydrosulfite oxidation. 



Table 13 

Regression Summary for Dependent Variable: ORP 
R= .99514836 R 2 = .99032025 Adjusted R 2 - .98709366 
F(i ,3)=306.93 p<.00041 Std.Error of estimate: 6.0576 



Intercept 



Ratio 
1460/1977 



BETA 



-.995148 



St. Err. 



of BETA 



.056803 



ORP=-206.467-75.062*I460/I977 



B 



-206.467 



-75.062 



St. Err. 



ofB 



t(3) 



9.147723 



-22.5703 



4.284508 



-17.5193 



p-level 



.000190 



.000405 



Table 14 

Correlations Hydrosulfite Oxidation 
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Marked correlations are significant at p < .05000 



TIME 



time 



1.00 



OR 
P 



.94 



1228 



1.00 



1460 



-.99 



1583 



1.00 



1977 



.94 



1998 



1.00 



I 

102 
2 



-.79 



Rati 
o 

583 
977 



-.94 



Rati 
o 

977 
998 



.90 



Rati 
o 

460 
977 



-.94 



Rati 
o 

977 
102 
2 



.90 



ORP 



1228 



1460 



1583 



1977 



1998 



11022 



Ratio 

583 

977 



Ratio 

977 

998 



Ratio 

460 

977 



Ratio 

977 

1022 



.94 



1.00 



-.98 



-.99 



-.95 



1.00 



-.99 



1.00 



.94 



1.00 



-.79 



-.94 



.90 



-.94 



-.98 



-.99 



1.00 



1.00 



-.95 



.97 



-.88 



-.96 



-.99 



.93 



1.00 



.90 



.96 



1.00 



-.94 



1.00 



1.00 



1.00 



1.00 



-.94 



.99 



.79 



.94 



-.91 



.94 



-.90 



.99 



.79 



.96 



1.00 



1.00 



-.94 



1.00 



.79 



.96 
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.96 
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-.96 
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.94 
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.94 



-.94 



.99 



.79 



.96 



-.89 



.96 
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1.00 



.79 
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-.89 



.95 



1.00 



-.95 



-.93 



-.87 



.98 



.95 



-.89 



-.87 



.99 



-.95 



-.93 



-.87 



.98 



-.87 



1.00 



.81 



.95 



-.91 



.94 



.81 



.95 



-.91 



.94 



-.91 



1.00 



.75 



-.91 



.76 



.75 



1.00 



-.77 



1.00 



-.91 



.76 



-.96 



-.77 



1.00 



-.81 



1.00 



-.77 



.99 



-.77 



1.00 



-.80 



.96 



-.90 



-.89 



-.89 



.99 



-.91 



-.96 



-.81 



.99 



-.80 



1.00 



The Oxidation/Reduction Potential (ORP) is a measurement of the potential for a reaction 
to occur. Oxidation-Reduction represents electron concentration and activity level. An 
ORP in the plus range indicates oxidation, i.e. the absence of energy, and an inability to 
perform additional chemical reactions. An ORP in the negative range indicates chemical 
reduction, i.e. the presence of electrons, potential energy, and the ability to generate 
additional chemical reactions. ORP is therefore a measure of energy potential. The more 
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negative the ORP, the more electrons present (in relation to the number of protons), and 
the more energy available. Biological redox reactions are a result of hydrogen being the 
essential electron donor, and oxygen being the essential electron acceptor. 

5 Oxidation reduction potential, ORP, measurements are used to determine the oxidizing or 
reducing properties of a solution. Typical applications include the neutralization of waste 
water containing sulfides, cyanides, chromates, nitrites or organic waste and controlling 
the addition of oxidants to drinking water, swimming pools or cooling towers. ORP can 
also be used to determine the ion activity of metals in solution and determine the endpoint 

1 0 of titrations. 

The standard method for determining the ORP is given in "Standard Methods for the 
Examination of Water and Waste Water", 18 ,h Edition, 1992, method 2580 B. Discussion 
of relevant background information is included in this reference. 



= 15 



The term potential of an oxidative reductive process within this specification is defined as 
the oxidation reduction potential and also as an oxidation reduction indicator. The 
potential of an oxidative reductive process is a measure of a composite state of the 
oxidation reduction capacity of a sample and not of the oxidative reductive capacity of 
20 the individual species. 

ORP is commonly measured as pE using a metal electrode, a reference electrode and a 
high input impedance millivolt meter such as a pH meter. However, electrochemical 
methods for measuring the ORP cannot discriminate between different 

25 oxidative/reductive species. In the environment or in industrial applications usually 
many different influences on the ORP from other oxidative or reductive species are 
present. In accordance with an embodiment of the present invention Raman spectroscopy 
is employed to determine the contribution to ORP of certain species more accurately by 
monitoring the oxidative/reductive properties and relative concentrations of the species of 

30 interest. Alternatively, additional oxidative/reductive species contributing to the ORP are 



39 



51-6 USCIP 



Patent 



monitored if their molecular vibrations are Raman active. 

The Nernst equation relates the oxidative and reductive properties of a certain regime. 
The Nernst equation can be expressed in pE notations as: 

pE = pE - + i logm Mod. [i] 

n {Area} 

pE is a notation related to the oxidation form and the reduction form of the species in a 
regime. Presented below is the example for a water regime. 

The limit of pE in water can be determined using the Nernst equation. The stability 
regime of water is determined by the reduction and oxidation of water. The reduction of 
water defines the lower limit to pE 

H + +e' ^-H 2 0 P] 
2 



an [3] 



pE^pE' +log 10 i 

fn2 



As pE' is zero by definition and the boundary condition commonly used is a hydrogen 
fugacity, f H2 of unity, this gives 

pE = -pH M 

The upper boundary for the stability of water is given by the oxidation reaction: 

pE = P E-+\og ]0 f ( %{H + } M 
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The boundary condition used is an oxygen fugacity of unity. The upper pE limit becomes: 



p 

1=; 

4* 



U 

5 

ru 

'bus? 



10 



15 



20 



pE = 2015 -pH [7] 

The pE for natural waters can be calculated using equation [6]. The unknowns which 
must be measured are pH and the 02 partial pressure (concentration). For example, the 
pE for surface water in equilibrium with the atmosphere (p G2 = 0.21 atm) and have a pH 
of 8 would be: 



pE = pE--pH + log l0 f$ PI 



^ = 20.75-8 + log 10 (0.21) X =12.58 [10] 

The stability boundaries for water were calculated earlier. The upper and lower limits, 
respectively, are: 

pE = 20.75 -pH t 11 ] 
P E = -pH H2] 

Measurements of Properties of Species Forming Complexes, Polymerized and 
25 Networked Structures, and Multiple degrees of Ionization. 

Figure 20 presents a series of silicate Raman spectra as a function of varying hydroxide 
concentration as taken from Prabir K. Dutta and Dah-Chung Shieh, published in Applied 
Spectroscopy, Vol. 39, No. 2, pp. 343-346 (1985). These spectra yield the intensities and 
ratios presented in Figures 21 and 22. Silicates can take different structural forms 
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through ionization and polymerization in dependence upon a pH value. Silicates can be 
used in a plurality of industrial applications. For example, silicates are added to a 
bleaching solution in order to stabilize peroxides. 

The set of spectra shows how silicate types vary with solution conditions. 

Figure 21 shows the intensities of the different peaks as derived from the spectra 
presented in Figure 20. Figure 21 shows a matrix plot of Raman intensities as a function 
of the HOVSi ratio as derived from Prabir K. Dutta and Dan-Chung Shieh, Applied 
Spectroscopy, Vol. 39, No. 2, 343-346 (1985). Peaks from silicate relate to different 
vibrational modes that depend on the ionization and degree of polymerization of the 
silicate. Vibrations from silicate monomers yield peaks at 925 cm" 1 (Si-O" stretch, 
monomer ionized) 772 cm*' (Si-O-H stretch, monomer not ionized), 482 cm* 1 (Si-OH 
stretch) and 446 cm' 1 (Si0 2 (OH) 2 2 ' symmetric bend). Dimer groups yield peaks at 598 
cm' 1 (Si-O-Si stretch, dimer bridge) and 1014 cm" 1 (Si03 stretch, dimer endgroup). 
Cyclic trimers have a breathing vibration at 53 1 cm" 1 . The peak at 1014 cm' 1 (Si03 
stretch, dimer endgroup) shifts to -1030 cm* 1 with cyclic trimers. 

The prior art has looked only at the Raman intensities. In accordance with the present 
invention ratios are used as they are capable of presenting different non-linear 
relationships with respect to solution properties than the intensities do. 

Ratios present the physical properties of the silicates better than intensities alone do. This 
is demonstrated in Figure 22 showing a matrix plot of Raman intensity ratios as a 
function of the HOVSi ratio from Prabir K. Dutta and Dah-Chung Shieh, Applied 
Spectroscopy, Vol. 39, No. 2, 343-346 (1985). It is noted that this article does not teach 
the use of ratios. The intensities given in this article were used to derive Raman intensity 
ratios in order to demonstrate that the use of Raman intensity ratios is much better in , 
predicting the physical properties of silicates. The ratios relate the relative concentrations 
of different functional groups characterizing the silicate speciation. These ratios provide 
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new variables that logically relate to the properties of the silicate solution. The ratio 
R579/1014 is related to the degree of polmerization or the chain length because 
individual intensities represent the amount of bridging/the amount of end groups, i.e. ratio 
of the peak intensity at 597 cm" 1 to the intensity at 1014 cm"'. The degree of 
polymerization decreases with increasing alkali concentration. The ratio R925/772 is 
related to the degree of amount of ionization of the silicate monomers (ratio of the peak 
intensity at 925 cm" 1 to the intensity at 772 cm' 1 ). The ratio R53 1/597 is related to the 
ratio of cyclic trimer groups to dimer groups (ratio of the peak intensity at 531 cm' 1 to the 
intensity at 597 cm' 1 ). The ratio R53 1/772 is related to the ratio of cyclic trimer groups to 
protonated monomer groups (ratio of the peak intensity at 53 1 cm' 1 to the intensity at 772 
cm' 1 ). The ratio R53 1/925 is related to the ratio of cyclic trimer groups to ionized 
monomer groups (ratio of the peak intensity at 53 1 cm 1 to the intensity at 925 cm 1 ). The 
ratio RIO 14/925 is related to the ratio of ionized dimer groups to ionized monomer 
groups (ratio of the peak intensity at 1014 cm' 1 to the intensity at 925 cm' 1 ). 

FT-Raman spectra were recorded using a Bruker IFS-88 Fourier Transform Infrared (FT- 
IR) spectrometer equipped with FRA 106 Raman accessory. Excitation was provided by 
an Nd:YAG laser (X - 1.064mm, v -9394 cm' 1 ). Solid state samples were measured 
using an aluminum sample holder and solutions were contained in a 0.5 cm cuvettes. 
Sepctra were recorded at a resolution of 4 cm' 1 , utilizing a laser power typical 50mw for a 
solid sample and 280mw for a solution sample. 

The spectrum in solid-state is different from the spectrum from the solution. In most 
cases an aqueous solution was measured. In general, Raman spectroscopy has an 
advantage over IR spectroscopy with respect to working on an aqueous phase. The 
detection limit can be pushed down to -0.05% at current experimental conditions (Laser 
power 280mw and single path at 180° sample collection configurations). In order to 
increase the detection limit, higher laser power and/or multi-path sample collection 
configuration are needed in FT-Raman spectroscopy. 
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Example of The Use of Raman Spectra and Ratios for Measurement and Control of 
Pulp Bleaching. 

There are 1 0 different samples (PH = 9.5,9.7, 10.0, 10.3, 10.5, 10.7 11.5, 11.7, 12.1 and 
12.7). Since each sample adjust the final pH by using buffer solution (pH = 6.9) and 1M 
NaOH, these 10 sample have been divided into two groups. There is one group with 
three peaks at 877, 990 and 1077 cm' 1 , and another group having two broad bands at 1 100 
and 1630 cm' 1 . These two broad bands are due to water and the Raman sample cell. As 
described above, water was subtracted from the obtained Raman spectra before 
processing the data. 

According to the reference Raman spectra of H 2 0 2 , the peak at 877 cm' 1 in bleaching 
water is due to the H 2 0 2 . The amount of the H 2 0 2 is approximately 0.1%. The peak at 
990 cm" 1 is due to S0 4 stretching band of the sulfate (reference see previous work). The 
amount of H 2 0 2 decreases with the pH increase, while the amount of the sulfate increases 
with the pH increase. The detection limit for H 2 0 2 and S0 4 can be lowered down to 0. 1% 
at current instrument condition. Interfering peaks due to carboxylic acid vibrations occur 
at 890 cm" 1 for acetic acid and 927 cm" 1 for the acetate ion. 

The third peak in bleaching water sample is at 1077 cm 1 . The intensity of this peak does 
not change with pH increase. 

Other Examples of Raman Spectra of Ultrafiltered Pulp Processing Waters 

Raman spectra of two white water samples from Millar Western Pulp (MWP) and one 
from Alberta Newsprint (ANC) are collected and shown in Figure 23. These samples are 
filtered through a 0.45 micron filter before data is collected. Only one band at 990 cm" 1 
appears in ANC water. It is the contribution of the S0 4 (sulfate ion) band from the 
sulfate. There are also some sulfates in two kinds of Millar Western Pulp water. The 
amount of sulfate is different according to the band intensities. The peak at 877 cm" 1 in 
both Millar Western Pulp water samples also indicates the existence of H 2 0 2 . Another 
band at 925 cm" 1 in two MWP water samples is possibly the contribution of the acetate 
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buffer solution. Figure 24 gives the Raman and IR spectra of acetic acid and acetic buffer 
solution (0.05M). The amounts of the buffer solution are different in these two kinds of 

different water samples. 

Other relevant small molecules and complex ions in the pulp and paper industry may be 
detected. These species include S0 4 > SO * H 2 0, C102, HC.03, silicates, acet.c acd, 
Chloric Acid HC103, Chlorate C103(-), Chlorous Acid HC102, Chlorite C102(-), 
Hypochlorous Acid HCIO. Hypochlorite CIO(-), phosphate, nitrate, nitrites 

Assignments from organic model compounds 

1 800-1 600 cm" region: The stretching bands of C=0 from -COOH groups were shown 
up in spectra of Glucuronic acid and Polygalacturonic acid. 

1500-1000 cm" 1 region: The CH bending bands of the relevant R groups are in this region 
except for one band at -1330 cm"' is due to the -OH bending. The CO stretching bands 
from different groups appear in 1200-1000 cm-1 region. 

1000-300 cm- 1 region: The bands appear from 1000-700 cm' 1 are due to different CH 
bonds wagging or scissoring modes. Most bands below 700 cm" 1 are due to molecular 
skeletal deformations. 

Raman spectrum of Vanillin in solid-state show two more CH stretching bands at 3077 
and 3030 cm" 1 belong to aromatic C-H bonds. Three bands at 1 662, 1 589 and 1 508 cm" 1 
are due to the different C=C bond stretching modes. These are main different 
characteristic bands between lignin and carbohydrates. PAS-IR spectrum of Kraft lignin 
is shown in Figure 2b. The strong -OH absorbent band and other hydrogen bonding 
absorbent bands in IR spectrum have strong negative effect on spectral quality. In order 
to avoid the disturbance from the -OH groups and hydrogen bonds in IR spectra, Raman 
data would be a good choice to work on. 
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Raman spectra of tristerain and abietic acid were obtained. These two compounds belong 
to the wood extractive categories. Raman spectra of these two compounds are typical fatty 
acids with different chain length. By simply comparing the intensities of v c . 0 bands at 
-1700 cm" 1 with the intensities of the CH stretching bands at 2800-3100 cm ', the chain 
length of abietic acid is shorter than tristerarin. 

The invention provides an additional set of process parameters that are closely tied to the 
process chemistry. Using these new variables with functions that predict, model and control 
the state and outcome of a multi-dimensional process provides a significant advantage over 
standard chemometric and linear regression methods. Advanced control systems including 
those utilizing fuzzy logic, time correlation analysis, neural networks, adaptive control, 
principle component analysis and partial least squares provides a means of developing 
software programming and logic solutions to problems with a high dimensionality. It is an 
object of this invention to provide new variables that may be used to better control input and 
output parameters in a complicated process with many interrelated variables. 

The above-described embodiments of the invention are intended to be examples of the 
present invention and alterations and modifications may be affected thereto, by those of 
skill in the art, without departing from the scope of the invention which is defined solely 
by the claims appended hereto. 
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